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Abstract 

After a brief review of the historical role of analyticity in the study of criti- 
cal phenomena, an account is given of recent discoveries of discretely holomorphic 
observables in critical two-dimensional lattice models. These are objects whose cor- 
relation functions satisfy a discrete version of the Cauchy-Riemann relations. Their 
existence appears to have a deep relation with the integrability of the model, and 
they are presumably the lattice versions of the truly holomorphic observables ap- 
pearing in the conformal field theory (CFT) describing the continuum limit. This 
hypothesis sheds light on the connection between CFT and integrability, and, if 
verified, can also be used to prove that the scaling limit of certain discrete curves 
in these models is described by Schramm-Loewner evolution (SLE). 

1 Introduction 

Analyticity has played a key role in the development of our mathematical understanding 
of the nature of critical phenomena. At various levels, the assertion of the analytic 
dependence of a suitable physical observable on some variable has provided a powerful 
starting point for theories, and its ultimate failure has informed the next step in their 
refinement. 

The first example of this was perhaps Landau theory, in which the free energy F[M] is 
postulated to be analytic in the order parameter M: 



F[M]= J ((VM) 2 + Tr(H • M) + r TrM 2 + A 3 TrM 3 + • • ■) d 



d x. 
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In this expression, all terms allowed by symmetry are in principal present. The require- 
ment that M should minimize F, 5F/5~M. = 0, implies that there are critical points (in 
the mathematical sense) at certain values of the parameters. The singular behavior close 
to these is described by the mathematical theory of bifurcations. The critical exponents 
which emerge from this analysis are super-universal in that they do not depend on the 
spatial dimension d or the symmetry of the order parameter. 

Of course, the successes and the limitations of Landau theory are well understood. It does 
not account for the fluctuations in the order parameter. In high enough dimensions these 
can usually be ignored as far as the universal behavior is concerned, but in general they 
should be taken into account by computing the full Landau- Ginzburg- Wilson partition 
function 

Z[B] = [ [dM(x)} e~ F ^ , 

J|q|<A 

with a suitable cut-off A on the short wavelength modes. This of course is impossible ex- 
cept in some simple cases, but these can often be used as a starting point for a perturbative 
expansion. However this approach generally fails near the critical point. 
The renormalization group (RG) of Wilson, Fisher and others was a way of dealing with 
this problem, by considering the response to changes in the cut-off A — > Ae~ e in such 
a way as to move the parameters into a region where perturbation theory is applicable. 
Requiring that the long-distance physics remains the same then leads to the RG flow 
equations 

f = - A t = -*«A» 

for the parameters {Xj} in F[M]. The critical behavior is then controlled by the fixed 
points where Pj({X}) = 0. One of the crucial assumptions of the RG theory is that these 
functions themselves are analytic in the {Xj} near the fixed point. Thus analyticity is 
moved to a higher level of abstraction: the non-analyticity of the free energy results from 
a potentially infinite number of applications of the analytic RG mapping. [Interestingly 
enough there is very little direct evidence for the strict validity of this assumption. It is 
largely based on analysis of perturbation theory which begs the question. Solvable models 
in two dimensions have so much analytic structure built into them that analyticity of the 
RG flows is almost inevitable, and possible detailed violations of analyticity have not, 
to my knowledge, been carefully studied in numerical simulations in higher dimensions. 
Indeed in one case, d = 0, when the beta-function can be evaluated analytically, it displays 
an essential singularity at the fixed point.] 

The next example of the use of analyticity is in the solution of integrable models in two 
dimensions. In this case the independent variable parameterizes the solution manifold on 
which the Boltzmann weights satisfy the Yang-Baxter equations. Its physical interpreta- 
tion, to be discussed later, is the degree of anisotropy of the model in two-dimensional 
euclidean space. The Yang-Baxter equations imply that the row-to-row transfer matrices 
at different values of this parameter commute. Assuming that the analytic properties of 
the local weights in the parameter lift to thermodynamic quantities, Baxter and others 
have shown that these obey functional relations which often determine them completely. 
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The other main approach to understanding two-dimensional critical behaviour attempts 
to describe the continuum critical scaling limit directly, and also exploits analyticity. This 
is the approach of conformal field theory (CFT), more recently also linked to Schramm- 
Loewner evolution (SLE) p]. In this case the analyticity is in the two-dimensional coordi- 
nates z = x + iy and z = x — iy. In CFT, correlation functions of certain observables are 
holomorphic functions of z (or antiholomorphic functions of z.) These observables are of 
two types: conserved currents corresponding to continuous symmetries, such as the stress 
tensor (T(z),T(z)), or so-called parafermions ip s (z) whose 2-point correlation function in 
the full plane have the form 

{if> 8 (zi)if>,(z2)) ~ (z 2 -z 2 )~ 2s , 
where s, the conformal spin, is in general fractional. 

These holomorphic (and antiholomorphic) objects are the building blocks of many CFTs. 
They also give an explicit meaning to the statement of conformal covariance: if we con- 
sider the scaling limit of a critical lattice model in a domain T>, then the behavior of a 
multi-point holomorphic correlator is completely determined by its singular behavior at 
coincident points and the boundary conditions on dT>. If these are themselves conformally 
covariant {e.g. if \im z ^QD aigip s {z) is determined by the boundary tangent angle), then, 
under any conformal mapping / : T> — > T>' of the interior of D to another domain V 

It is the purpose of this paper to review recent work which forges a link between these last 
two realizations of analyticity More specifically, starting from certain lattice models we 
identify so-called discretely holomorphic observables, whose correlators satisfy a lattice 
version of the Cauchy-Riemann equations. These have fractional spin by construction, 
and are presumably the lattice precursors of the parafermions in the corresponding CFT. 
We find, as expected, that discrete holomorphicity holds only when the Boltzmann weights 
lie on the critical manifold of the model, but, more surprisingly, that they also lie on the 
integrable critical manifold, that is, they satisfy the Yang-Baxter equations. 
This work has been described in detail in several papers. In [21 [3] a discretely holomorphic 
observable was identified for the random cluster representation of the g-state Potts model. 
In [1] several were found, directly from the Boltzmann weights, for the clock models, 
and the close connection to integrability was first observed. Smirnov [3] also found such 
an observable for the random curve representation of the O(n) model on the honeycomb 
lattice. This was extended to the O(n) model, and further generalizations, on the square 
lattice in [5], where the parameter space is large enough to make clear the connection to 
integrability. Smirnov's objective in identifying these holomorphic objects associated with 
curves was to show that their scaling limit in described by Schramm-Loewner evolution 
(SLE). This approach is fully described in [3] and we refer the reader there. Suffice it 
to say that proving that the discretely holomorphic lattice observables go over into fully 
holomorphic quantities in the scaling limit is in general a very difficult problem which has 
been solved completely only in a few special cases, in particular that of the Ising model 
(q = 2 orn = 1) [3]. 
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Figure 1: A face of the square lattice. The discretely holomorphic function F(z) is defined 
at the mid-points of the edges. 



In this account we shall describe in detail only two examples. The first is the simplest: 
the case N = 2 of the Zn model, more usually known as the Ising model. In this case 
the parafermion has conformal spin s = \ and is related to the fermionic objects used 
in the various exact solutions of the model, going back to Onsager. It is defined, in a 
standard way, as a product of nearby order and disorder variables. This identification 
dates back (at least) to Fradkin and Kadanoff [6], but the simple argument given in [I] 
that a particular lattice definition is discretely holomorphic is, to our knowledge, new. An 
essential ingredient in the definition of a disorder operator is that the model must possess 
a dual description in the sense of Kramers and Wannier. However, the second example is 
the 0(n) model which does not possess this. In this case the parafermionic observables 
are instead defined in terms of the random curve representation of the model, which is 
nonlocal in terms of the original spin degrees of freedom. 

Before proceeding, however, we give a proper definition of discrete holomorphicity. Sup- 
pose Q is a planar graph embedded in R 2 , for example a square lattice. Let F(zij) be 
a complex- valued function defined at the midpoints of each edge (ij). Then F is 
discretely holomorphic on Q if 

£ F(z lJ )(z J -z i ) = 0, (1) 

where the sum is over the edges of each face T of Q. This is a discrete version of the 
contour integral. For a square lattice (see Fig. [1]) it reduces to 

F(z 12 ) + iF(z 23 ) + i 2 F{z u ) + * 3 F(24i) = . 

A little thought shows, however, that the total number of such equations, one for each 
face of Q, is in general far less than the number of unknowns, one for each edge. Therefore 
even on the lattice this system does not determine F(zij) for suitable boundary conditions, 
unless further information is available. (For the Ising model it turns out that more is 
known about the phase of F, so the system is rigid.) However, in the continuum limit 
as the lattice spacing tends to zero, we can approximate arbitrarily closely the contour 
integral J c F(z)dz by a sum of the left-hand side of ([T|) over faces T which tile the interior 
of C, and therefore assert that is vanishes for all reasonable contours C. Morera's theorem 
then assures us that if F(z) is continuous then it is analytic. In field theory the continuity 
of correlation functions is usually taken for granted, but strictly this needs to be proved. 
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Figure 2: The string defining a disorder operator. 




s(rj 



\y-rn 



\i(R 3 ) 



s(rj) 



Figure 3: Elementary square on whose edges the parafermionic observables are defined. 

2 The Ising and Zn models 

The square lattice Ising model has spins s(r) = ±1 at the vertices r. The Boltzmann 
weights are exp(— H) with H = —J2 r r'Jrr' s ( r ) s ( r> ), where the sum is over edges (rr'). 
Note that the weights can also be written 

Y\ rr i (1 + (tanh J rr i)s(r)s(r')). The insertion of a disorder operator n(R) at the dual vertex 
R corresponds to changing J rr i — » —J rr i on all edges which cross a 'string' attaching R to 
the boundary, see Fig. [2j That is 



The correlator (fi(Ri) fi(R 2 )) corresponds to such a string connecting Ri and R 2 , and is 
invariant under deformations of the string. 

The parafermionic variables ip s {rR) are defined on the midpoints of each edge (rR) con- 
necting the vertex r to a neighboring dual vertex R, which form the covering lattice: 



n 



1 — (tanh J rr i)s{r)s(r') 



(2) 



1 + (tanh J rr t)s{r)s{r') 



ijj s (rR) = s(r) ■ fi(R) e 



-is9(rR) 



Here 8(rR) is the angle that (rR) makes with (say) the positive x-axis, but we have to 

be careful about its multivaluedness (see below). 

Consider now an elementary square (see Fig. [3]). From ([2]) we have 



(1 + (tanh J y )s( ri )s(r 2 )) fi(R 4 ) = (1 - (tanh J y )s( ri )s(r 2 )) fi(R 3 ) , (3) 

where we have set J rr < = J XtV depending on whether (rr') is parallel to the x or y-axis. 
Now multiply ([3]) by s(ri) and s{r2) and use s(r) 2 = 1. This leads to two linear equations 
in the four parafermionic variables defined on the edges of the square. Simple algebra 
then shows that these imply the discrete holomorphicity condition (CQ) as long as: 

• we distort each square into a rhombus whose angle depends on the anisotropy J y j J x ; 

• we are careful, since s is in general non-integer, to define 9{rR) consistently so that 
it varies only by increments in the interval (— tt, it) on going around the square; 

• the couplings lie on the critical manifold sinh J x sinh J y — 1 . 

This example illustrates the general results discussed in the introduction. It also shows 
how the lattice should be embedded in R 2 so that its continuum limit is rotationally 
and conformally invariant. However since the nearest neighbor Ising model is integrable 
for all values of the couplings (J x , J y ), it is not general enough to illustrate the role of 
integrability. 

This is afforded by the generalization to the models [I]. In these models the Ising 
spins are generalized to complex roots of unity such that s(r) N = 1. The Boltzmann 
weights take the form 

/ AT-l \ 

n i+E^MMO) fc +c.c. . 

rr' \ k=l / 

It is then found that one can identify discretely holomorphic parafermions, again the 
product of neighboring order and disorder operators, with spins 

, = M^ (l<*<*/2), 
as long as the parameters lie on the manifold 

{k) _ ^ sm((7rj + a)/N) 



x x .' [ a 



A i sin ((tt(; + 1) - a)/N) [ ] 



and x y k \a) = x y k '(7r/2 — a). Here a is the half-angle at the vertex of the rhombus 
into which the square lattice must be distorted, and therefore measures the degree of 
anisotropy. However this manifold is precisely the critical integrable case for the general 
nearest neighbor model found by Fateev and Zamolodchikov [7j. The values of the 
conformal spins above agree precisely with those of the parafermionic holomorphic con- 
formal fields in the CFT postulated by the same authors [8] to describe the continuum 
limit of this model. 

The above analysis can be generalized simply to a general graph £ whose faces are 2- 
colorable, sometimes called a Baxter lattice - see Fig. HI The Z N spins s(r) are defined 
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Figure 4: Part of a Baxter lattice. The faces of the graph C formed by the curved lines 
are 2-colorable (not shown). Order variables s and disorder variables fi are associated 
with alternately colored faces respectively. The covering lattice is shown as solid lines. 
The theorem of Kenyon and Schlenker [9] asserts that for every such graph the covering 
lattice admits a rhombic embedding in the plane, that is one where all its edges have the 
same length. 



on the faces of a given color, and the disorder variables fi(R) on the others. Neighboring 
order and disorder operators then lie at the vertices of quadrilaterals which tile the plane. 
A theorem due to Kenyon and Schlenker [5] asserts that, under rather general conditions, 
they can can be distorted in the plane so they are all rhombi, that is, they form an 
isoradial lattice, on which nearest neighbors are all the same distance apart. If now the 
interactions across each rhombus are chosen to satisfy (jl]), where a is the half-angle of 
the rhombus, then the associated parafermion is discretely holomorphic on the isoradial 
lattice. Moreover the corresponding weights satisfy the star-triangle, or Yang-Baxter 
equations, as explained in Fig. 

We note that although we have defined the parafermions in this case directly in terms of 
the modified Boltzmann weights, they may also be identified with observables of random 
curves in the model. For the case of the Ising model it is well known that the high- 
temperature expansion, in powers of tanh J x y , of the correlation function (s(r)s(O)) can 
be expressed as sum over graphs on the lattice. The argument is clearer for the honeycomb 
lattice, when these graphs consist of non-intersecting closed loops and one open curve 7 
from to r. The additional complications on the square lattice are believed to be irrelevant 
to the scaling limit. The introduction of the disorder operator yu(-R) neighboring s(r) acts 
to weight configurations of the open curve with (— 1)^ where N is the number of times 7 
crosses the string. This can also be written as e~ ls9 ° r where 8 0r is the winding angle of 7. 
In the scaling limit this takes arbitrarily large values and so it does not matter whether 
we compute it in increments of ±2ir as it crosses the string, or simply in increments of 
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Figure 5: Two different tilings of a hexagon by the same set of three rhombi. The right- 
hand case has an additional vertex associated to an order operator s as compared to that 
on the right. Discrete holomorphicity for each rhombus fixes the couplings on the dashed 
lines to be related by the star-triangle transformation. The two pictures are also related 
in the original graph £ by moving one of the curves past the vertex formed by the other 
two - the Yang-Baxter relation. 



±| as we proceed along 7. This latter definition corresponds to the one used by Smirnov 

3 The O(n) model 

As a second example we consider the so-called 0(n) model on the square lattice, first 
considered by Nienhuis [10J. Although it can be written in terms of n-component spins 
s a ( r ) (a = 1)2, ... ,n) located on the edges of a square lattice through local, but non- 
nearest neighbor, interactions, it is more easily formulated in terms of a gas of dilute 
non-intersecting planar loops. Each elementary square can take one of the configurations 
shown in Fig. El with the indicated weights. When patched together these curves form 
closed loops, each of which receives an additional weight n, which no longer has to be 
an integer and may be parameterized by n = — 2 cos 27/ with < rj < |. Note that 
we have grouped the anisotropic weights so that there is symmetry under reflections in 
the diagonal axes. This symmetry is preserved when the plaquettes are deformed into 
rhombi. Since every loop configuration has an even number of plaquettes of type u\ or 
it2, the change (u\, U2) — > (—Mi, — M2) does not affect the Boltzmann weights. 
In order to define a suitable candidate for the parafermionic observable we must consider 
not only closed loops but also open curves which begin and end at different points, say 
and r. Such configurations arise, as in the previous section, if we compute the correlator 
(s a (r)s a (Q)} of the O(n) spins. Once again, to define an object with non-zero spin, we 
additionally weight each open curve by a phase factor e~ ls0Or , where #0 is the winding 
angle from to r, the accumulation of the turns through ±|. Note that this can be 
arbitrarily large. This type of observable for the O(n) model was first considered on the 
honeycomb lattice by Smirnov [3]. 
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Figure 6: Vertices of the 0(n) loop model on the square lattice. 



Now consider the contributions to the holomorphicity equation (TjQ) where r is one of the 
edges of an elementary square. Since the open curve ends at r, it must first intersect the 
square on some edge r' (not necessarily the same as r). The configurations can then be 
decomposed according to the different possible r'. Without loss of generality we take it 
to be on the lower edge in Fig. [7J The configurations can then be further decomposed 
into classes corresponding to the various cases shown in this figure. In the last three cases 
a curve connects two of the other edges: it may be part of a closed loop or of the open 
curve. The idea is to sum over all possible ways of connecting up the curves internally 
through the square, keeping the external configuration fixed, and to try to satisfy ([T|) for 
each of these configurations. This can be done, and yields the following linear system for 
the weights: 



t + fiUi — a\~ u 2 — v 
-X^ui + nu 2 + A/if — /iA~ 1 (wJi + nw 2 ) 
nu\ — \u 2 — /iA~ 2 t> + jj(nwi + w 2 ) 
-fi\~ 2 Ui + fi\u 2 + nv — \~ 2 wi — \ 2 w 2 










(5) 
(6) 
(7) 



where we have set A = e t7TS , ip = (s + 
complex linear equations for six real unknowns 
(t, u\, u 2 , v, wi, w 2 ), and we have the relations: 



e llf . For real weights, ([5H8]) are four 



Im 



(n + 1) {[5]} -A//- 1 {EH 



-1 



Im 



\fi-\X 2 - n\~ 2 ) © + fi-\n\ 2 - \- 2 ) {[3} 



in 




0. 



Thus, we can generally reduce ([5H8]) to a 6 x 6 real system. 

There are two classes of solutions, for vanishing and non-vanishing v. First, if v = 0, then 
the configurations corresponding to {EH never occur, and so this equation does not hold. 
In the special case n = 1, there exists a non-trival solution for any value of s: 



t = sin 7is , Ui = sin(<yj — ns) , u 2 = sin (p , Wi + w 2 



sm7rs 



(9) 



It turns out [5] that this model can be mapped onto the six- vertex model (see Figure El), 
with weights uj\ = uj 2 = sin(<£> — its), uj^ = 004 = sin (p,uj 5 = ljq = shirrs. The corresponding 
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Figure 7: Loop configurations with one edge of the face T connected to point 0. 
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Figure 8: Mapping of the 0(n) model onto the six-vertex model for n = l,v = 0. 



anisotropy parameter is A = cos ns. This is an example of a model admitting a holomor- 
phic observable on the lattice, but for which the scaling limit of the corresponding curve 
cannot be described by simple SLE. This is because, for ordinary SLE, the central charge 
of the CFT is directly related to the SLE parameter k [1] and hence to the conformal 
spin s: c = 2s(5 — 8s) /(2s + 1). In the present case, since the boundary conditions for 
the six- vertex model are not twisted, its scaling limit has central charge c = 1 for all A. 
However the conformal spin s varies continuously with A. Therefore the scaling limit 
of the curve can be SLE, with k — 4, for at most one value (in fact A = l/y/2.) We 
conjecture that other values of A in fact correspond to a variant of SLE called SLE(4, p). 
For v — , n ^ — 1, we get a 5 x 5 linear system, with determinant (n 2 — l) 2 sin <p sin((p — 
its). Imposing sin<^ = yields (s + l)a = ran and in turn s = m', where m,m' are 
integers. Thus, for the solution to exist at any value of a, we have to set s — — 1. The 
Boltzmann weights are then: 

t = -Mi - U 2 , W\ = -Mi , W 2 = —U 2 ■ (10) 

The solution of the case sin(y9 — ns) = is similar, and leads to the same Boltzmann 
weights and spin s = —1. If we change the sign of Ui, u 2 , then the model (fTUI) is equivalent 
to a dense loop model which corresponds to the critical 9- state Potts model, with weight 
per loop yfq = n + 1. To see this, fill empty spaces with loops of weight 1 (ghost loops). 



10 



The local weights do not depend on the type of loops involved (actual or ghost loops), so 
each loop has an overall weight n + 1. As a consequence, the dense loop model has a lattice 
antiholomorphic observable (s < 0), besides the holomorphic one found in [2j. However, 
several arguments rule out the hypothesis that this corresponds to an antiholomorphic 
field in the continuum limit. First, ip s =-i(z) is lattice antiholomorphic for any Q > 0, 
whereas it is well known that the self-dual Potts model is only critical for < Q < 4. 
Furthermore, the ratio Ui/u 2 in [10] does not depend on the angle a, which means that 
the same model has an antiholomorphic observable for any deformation angle: this is not 
acceptable physically in the continuum limit. So we conclude that, in the case v = and 
generic -1, the holomorphicity conditions (CQ) for the dilute O(n) model merely lead 
to the case of the dense loop model, but the corresponding ip s {z) is not a candidate for 
an antiholomorphic field in the continuum limit. 

Let us now discuss the solutions of second class (v ^ 0), for a generic value of n. We get 
the 6x6 real system: 

(Re © , Re © , Im © , Re {JZj}, Im ©, Re ® ) , 

with determinant: (n 2 — 1) sin<^sin((^ — its) (2 cos47rs — 3n + n 3 ). Non-trivial solutions 
exist if the spin satisfies: 

cos 47ts = cos 6r] . (11) 

The various solutions to [11] can be parameterized by extending the range of rj to [— tt, it], 
and setting: 

s = ^-~. (12) 
2vr 2 V ; 

Then, we get the second class of solutions, with Boltzmann weights: 

t = - sin (2</?- 3?7/2) + sin5r//2-sin3r//2 + sinr//2 (13) 

u\ = — 2 sin r/ cos (377/2 — <p) (14) 

U2 = —2 sin rj sin ip (15) 

v = — 2 simp cos (3r]/2 — ip) (16) 

wi = — 2 sin(p — rj) cos (3?]/2 — <p) (17) 

W2 = 2 cos (77/2 — <p) sinp . (18) 

A remarkable fact is that the weights f[T3T[T8l) are a solution of the Yang-Baxter equations 
for the O(n) loop model on the square lattice. Indeed, after a change of variables p — > 
ip + (it + 77) /4, they coincide with the integrable weights in [TDJ . So, by solving the 
holomorphicity equations (J3HS]) on a deformed lattice, we recover the integrable weights. 
Other more complicated loops models (for example one with different types of loops known 
as the C2(l) model [11]) can also be studied with similar results [5]. 

4 Conclusions and further remarks 

We have given two main examples of lattice models, the Zm model and the O(n) model 
on a square lattice, in which observables can be identified whose correlators are discretely 
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holomorphic, as long as the weights are both critical and satisfy the Yang-Baxter relations. 
This is surprising, since the holomorphicity conditions are linear in the weights, and work 
for a fixed value of the anisotropy parameter, while the Yang-Baxter relations are cubic 
functional equations for the weights. While in the case of the Zn models the connection 
to integrability may be understood by generalizing the problem to an inhomogeneous 
Baxter lattice, this explanation is at present missing for the loop models. While it would 
be nice to elevate these observations to a more general result connecting holomorphicity, 
integrability and conformal field theory, the counter-example given in Sec. [3] in which a 
lattice holomorphic observable apparently does not correspond to a conformal field should 
warn us that there may be subtleties. For the models for larger values of N, problems 
in a lattice identification of the value of the conformal spin s were also noted in [4]. 
Although it is to be hoped that the results of Smirnov [3] in using these holomorphic 
observables to prove that the scaling limit of lattice curves is given by SLE can be extended 
to other models, the examples we have given which correspond to CFTs with central 
charge c > 1 show that this may not always be the case. 

Finally, it is be hoped that at some Statistical Mechanics Conference in the future the 
correct extension of these ideas to higher dimensions will be announced! 
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